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Abstract 

The Lagrangian and Hamihonian structures for an ideal gauge-charged fluid 
are determined. Using a Kaluza-Klein point of view, the equations of motion are 
obtained by Lagrangian and Poisson reductions associated to the automorphism 
group of a principal bundle. As a consequence of the Lagrangian approach, a 
Kelvin-Noether theorem is obtained. The Hamiltonian formulation determines a 
non-canonical Poisson bracket associated to these equations. 



1 Introduction 



The equations of motion of an ideal incompressible fluid on an oriented Riemannian 
manifold (M, g) are given by the Euler equations 



dv _ , 
— + VyV = -gradp, 
ot 



where the Eulerian velocity f is a divergence free ve ctor field, p is the pressure, and V is 



the Levi-Civita covariant derivative associated to g. lArnoldl |1966| has shown that equa- 



tions (11. ip are formally the spatial representation of the geodesic spray on the volume- 



on M. See also 


Abraham and Marsden 


1978 . 


Ebin and Marsden 


1970 


give the analytic for- 



mulation and many rigorous results concerning the Euler and Navier-Stokes equations 
derived from this geometric point of view. From the Hamiltonian perspective, equations 
f ll.ip are the Lie-Poisson eq uations on the Lie algebra Xniy fAf) of ©^(M), consisting of 



divergence free vector fields (IMarsden and WeinsteinI |l983| ). Here, Xdiv(^) is identified 



with its dual by the weak L^-pairing. 
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In iMarsden et. al.l 1983l | and iMarsden. Ratiu. and WeinsteinI 1984| . this approach 
is generahzed to the case of the motion of an ideal compressible adiabatic fluid 



dv 
di 



+ V„f = - gradp, 



dp 

+ divipv) = 0, 



1.2) 



ds 
di 



+ ds(v) 



0, 



where p is the mass density, s is the specific entropy, and p is the pressure. In this 
case, the configuration space is the full diffeomorphism group T>{M) and equations (11.21) 
are obtained via Li e-Poisson reduction for seraidirec t products. The Euler-Poincare 
approach is given in lHolm. Marsden. and Ratiul 1998 . 

In this work we generalize the two previous procedures to the case of a classical 
charged ideal fluid. More precisely, using a Kaluza-Klein point of view, we obta in the 
equations of motion by Lagrange- Poincare fsee lCendra. Marsden. and Ratiul |2001| ) and 



Poisson reduction by a symmetry group (see, e.g. IMarsden and Ratiul |1994| . §10.5). We 
consider on M a G-principal bundle P ^ M and enlarge the configuration space from 
the group of diffeomorphisms of M to the product of the group of automorphisms of P 
with the field variables. 

U G = we recover the Euler-Maxwell equations describing the motion of an 
electrically charged fluid 



^ dv „ a ,^ 

+ V„t; = — E + f X B 
m 



1 



dt 

dp 

di 
ds 

di 
dE 

'dt 



■ gradp. 



+ div{pv) = 0, 
+ ds(v) = 0, 



1.3) 



curl B 



m 



-pv, 



dB 

'dt 



curl E, 



div E = — p, div B = 0, 

m 

where v is the Eulerian velocity, E is the electric field, B is the magnetic field, m is 
the mass of the charged fluid particles, and the constant q is the electric charge of the 
particles. The Hami ltonian structure of th e incompressible Euler-Maxwell equations is 



already presented in IMarsden et. al.l |1983 



Returning to the general case of a G-principal bundle, we will show that the Lagrange 
Poincare and Poisson reduction methods lead to the equations for an ideal compressible 



adiab atic fluid carrying a gauge-charge, as given in iGibbons. Holm, and Kupershmidt 



1983|. We call these equations the Euler- Yang-Mills equations. The physically rele- 
vant examples are obtained for G = SU(2) or G = SU(3) in which case the associated 
fluid motion goes also under the name of chromohydrodynamics. For a Lagrangian 
description of the Euler- Yang-Mills equations and the associated variational principle 
formulated in local coordinates both in the non-relativistic and relativistic versions. 
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see 



Bistrovic. Jackiw. Li. Nair. and Pil 2003 



2004 



and iJackiw. Nair. Pi. and Polychronakos 

In these papers the variations are constrained ac cording to the general Euler- 

Poinc are variational principle for field theories presented in lCastrillon-Lopez. Ratiu. and ShkoUer 
2000| . 



The physical interpretation of the equations obtained by the methods given in this 
paper is the following. The evolution of the fiuid particles as well as of the gauge-charge 
density of the fiuid is given by a curve ipt in the automorphism group of a principal 
bundle P — > M. In fact, ipt is the fiow of a time- dependent vector field Ut on P. 
This vector field induces a time-dependent vector field Vt on M, which represents the 
Eulerian velocity of the fluid. Given the evolution of the Yang-Mills fields poten- 
tial and of the mass density, the vector field Ut induces also a Lie algebra valued and 
time-dependent function which represents the gauge-charge density of the fluid. 
Note the analogy with the classical Kaluza-Klein construction appearing in the formu- 
lation of the equations of motion for a charged particle in a Yang-Mills field. See also 



Bambah. Mahajan. and Mukkul [2006l | who generalize the Eulerian fiuid velocity vector 
to include a non-Abelian, or gauge, index. 

The paper is organized as follows. To fix notations and conventions, we summarize 
in [|2] some basic facts about principal bundles, connections, automorphisms, and gauge 
groups. The Hamiltonian and Lagrangian formulations of the Maxwell equations are 
recalled in ^ and generalized to the case of the Yang-Mills fields equations. The La- 
grangian formulation of the motion of a charged classical particle in a Yang-Mills field, 
that is, the Wong equations, are presented in §11 In ^ it is shown that the compressible 
and incompressible Euler- Yang-Mills equations consist of coupled equations. These are 
the Euler-Poincare equations of a semidirect product (associated to the automorphism 
group of a principal bundle) for the fiuid and charge variables and the Yang-Mills equa- 
tions for the vector potential (that is, the connection) and the "electric part" of the 
Yang-Mills field. The Hamiltonian counterpart of this result is presented in ^ one 
obtains coupled equations consisting of Lie-Poisson equations on the same dual for the 
fiuid and charge variables together with the Yang-Mills equations. Formally the Gauss 
equation relating the gauge-charge and the "electric part" of the Yang-Mills field is 
missing from this system. However, it is obtained by conservation of the momentum 
map associated to the invariance under gauge transformations. We naturally obtain the 
non-canonical Poisson bracket associated to the Hamiltonian formulation of the Euler- 
Yang- Mills equations. By applying the general process of reduction by stages, we recover 
some already known results about the Eu ler-Maxwell equations. We also show tha t the 



tw o different Poisson bra ckets derived in iGibbons. Holm, and Kupershmidt 
in 



1983 and 



Marsden et. al 



19831 ] are in fact obtained by Poisson reduction, at different stages, 
of the same canonical Poisson structure. Finally, in ^ we present a Kelvin-Noether 
Theorem for the Euler- Yang- Mills equations. 
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2 Connections, automorphisms, and gauge transfor- 
mations 



In this section we recall basic notions related to principal bundles. We also introduce 
notations and conventions that will be used throughout the paper. 



2.1 Principal and adjoint bundles 

Consider a smooth free and proper right action 

of a Lie group G on a manifold P. Thus we get the principal bundle 

n-.P^M:^ P/G, 

where M is endowed with the unique manifold structure for which tt is a submersion. 

To any element ^ in the Lie algebra of G there corresponds a vector field on P, 
called the infinitesimal generator, defined by 

At any p & P, these vector fields generate the vertical subspace 

VpP := {^p{p)\^eQ}^ker{T,n). 

Recall that the adjoint vector bundle is 

AdP := P Xg0 ^ M, 

where the quotient is taken relative to the right action {g, (p, C)) ^ {^gip)^ -^^g-^iO)- 
The elements of AdP are denoted by [p, for (p, ^ P x Q- There is a Lie bracket 
operation [•, -J^ on each fiber (AdP)^ depending smoothly on a; e M; it is defined by 

[\p,^]g,\p,v]g\^ ■= \p,[^,v]]g 
for \p, C]g, [p, v]g e Ad P, 7r(p) = x. 

2.2 Exterior forms on adjoint bundles 

Consider the space Q,''{P, q) of g- valued /c-forms on P and let Q'=(P, q) be the subspace 
of Q'^(P, g) consisting of g- valued /c-forms u; such that: 

(1) ^*gUJ = Adg-lOUJ, 

(2) if one of Ui, Uk € TpP is vertical then uj{ui, ■■■,Uk) = 0. 
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The real vector space f2'^(P, g) is naturally isomorphic to f2^(M, AdP), the space 
of Ad P- valued /c-forms on M. Indeed, to each uj G VL^{P,q) corresponds a /c-form 
uj G f2^(M, Ad P) whose value on f i, .., t>fc G T^jM is given by 

Uj{x){vi,...,Vk) := ...,Mfc)]G, (2.1) 

where p G P is such that Ti{p) = x and Ui G TpP are such that Tp7c{ui) = Vi. 

To define the inverse of the map ~: VI^{P,q) Vl''{M,MP) we introduce first for 
every p E P the R-linear map ip : (AdP)^; Q, x := tt{p) G M, by 

(b> ii]g) ■= ^, where ^ is such that [p, = [q, vh; 

in this formula p,q E P and ^,f] E Q- Equivalently, this definition can be restated as 

ip {[q,ri]G) := Adg?7, where G G is uniquely determined by g = 

Then the definition of the Lie bracket on each fiber (AdP)^ of the adjoint bundle AdP 
immediately implies that ip : (AdP)^ — g is a Lie algebra isomorphism. In addition, 
= Adg-i oip for every p E P and g E G. Using the maps ip for every p E P, define 
the inverse of ~: lF(P,0) n''{M,AdP) by 

uj{p){ui,...,Uk) ■.= ip(u{TT{p))(Tpn{ui),...,Tp7i{ui))^, (2.2) 

for any p E P and Ui,...,Uk G TpP. The identity «$g(p) = Adg-i oip ensures that 

ijen^iP^. 

Since l)0(P,g) = -Fg(P,0) := {/ : P ^ g | / o $^ = Ad^-i of} and n%M,AdP) = 
r(AdP), the space of sections of AdP, we shall use the notations f2°(P, g) = J^g{P,q) 
and 11° (M, AdP) = r(AdP) interchangeably. We have hence J^g(^,0) = r(AdP) as 
Lie algebras, the isomorphism / G J-'g{P,q) i-^ / G r(AdP) being given by fl2.ip . that 
is, 7(7r(p)) = [pJ{p)]g- 

2.3 Connections and covariant differentials 

A principal connection on P is a g- valued 1-form A G ^^{P, g) such that 

^*gA = Adg-i and A{^p) = ^ 

The set of all connections will be denoted by Conn{P). It is an affine space with underly- 
ing vector space fi^(P, g). Recall that a connection induces a splitting TpP = VpP®HpP 
of the tangent space into the vertical and horizontal subspace defined by 

HpP := kei{A{j))). 

The covariant exterior differential associated to A is the map d"^ : ^'^(P, g) — »• 
fi'=+i(P,g) defined by 

d-^uj{p){ui, ...,Uk) ■■= duj{p) (horp(Mi), horp(Mfc)) , 
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where hoYp{ui) is the horizontal part of G TpP, i = 1, . . . , k. Note that for uj G fl''{P, q) 
we have d'^u G i7'=+^(P, g). For / G J^g{Pi fl) and u; G fl) we have the formulas 



d-^/(p)(n) = df{p){u) + /(p)] and 

d-^uj{p){u,v) = diu{p){u,v) + [A{p){u),uj{p){v)] - [A{p){v),uj{p){u)] 

for any u,v E TpP. 

The curvature of the connection A is, by definition, the 2-form 

B:=d-^AeW{P,Q). 

The curvature B verifies the Cartan Structure Equations and the Bianchi Identity given 
respectively by 

B{u,v) = dAiu,v) + [Aiu),Aiv)] and d-^-B = 0. 

The following lemma will be useful for future computations. 

Lemma 2.1 Let A G Conn{P), B its curvature, and f G J-'g{P,q)- Then 

d^d-^f{u,v) = [Biu,v),f]. 

Proof. For any U,V E X{P) we have 

d-^{d-^ f){U,V) = d{d^ f){U,V) + [A{U),d-^f{V)] - [A{V),d-^f{U)] 
= d (d/ + [A, /]) (U, V) + [AiU), dfiV) + [AiV), /]] 

-[A{V),df{U) + [A{U),f]] 
= + d{[A{V), mU) - d{[A{U), mv) - [A{[U, V]), f] 

+ [AiU), dfiV)] + [AiU), [AiV), /]] - [AiV), dfiU)] 

- [AiV), [AiU), f]] 
= [diAiV))iU),f] - [diAiU))iV),f] - [Ai[U,V]),f] 

+ [AiU),[AiV),f]]-[AiV),[AiU),f]] 
= [dAiU,V),f] + [[AiU),AiV)],f] 
= [BiU,V),f]. m 

Recall also that a principal connection ^ on P induces an affine connection and a co- 
variant derivative, denoted respectively by V-^ and on th e vector bundles Ad P — » M 
and (A dP)* —>■ M (see e.g. iKobayashi and Nomizul 1963| or lCendra. Marsden. and Ratiu 
200l| l 



Given a Riemannian metric on M and a connection A on P, we can define the 
covariant codifferential 



S^:nkiP,Q)-^nk-^iP,g); 



see, e.g.. Definition 4.2.8 in iBleeckerl [19811 
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2.4 Bundle metrics 



We assume throughout this paper that the Lie algebra g has a distinguished inner 
product 7 satisfying 

7(Adg ^, Adg T]) = 7(^, T]), for all g eG and all ^,r] ^ 9- 

Such an inner product is said to be Ad-invariant and satisfies the relation 

7([C,e],^) + 7(e,[C,^]) = for all ^,vX e Q. (2.3) 

For example, if G is compact such an inner product always exists. If G is reductive one 
can always find such a nondegenerate 7 but it may be indefinite. 

Given a Riemannian metric g on M and an Ad-invariant inner product 7 on g we can 
define a Riemannian metric (77 on the vector bundles A'^(M, Ad P) —>■ M of Ad P- valued 
exterior fc-forms on M. Indeed, the inner product 7 induces a Riemannian metric on the 
vector bundle AdP M whose value on ^Jc, [v^vIg ^ (AdP)i., x := 'nip), J9 G P, is 
given by 

7x([p,^]g, [p,^]g) := 7(^,^)- 

Denote, by abuse of notation, by the same letter 7 the smooth vector bundle metric on 
AdP defined by 7|(AdP)^ '■= Ix- Let 'g denote the Riemannian metric induced by g on 
the vector bundles A'^M — >■ M of exterior fc-forms on M. The Riemannian metric (77 
on the vector bundle A'^(M, AdP) ^ M is then constructed in the following manner. 
If I3x e A'^(M, Ad P)a;, write = a"' fa and = P'^fa, where {/„} is a basis of the 
fiber (AdP)^, and a^^p"" E {A.^M)^. Then define 

{9l)x{ax,Px) := 7a6^(a'',/5''), 

where 7^6 := ■jxifa, fb)- It is easy to verify that this construction is independent on the 
choice of the basis in each fiber (AdP)^;. 

Let M be a compact oriented boundaryless manif old. If a G r2 ^(M, AdP) and 



p e n''+^{M,AdP) we have (see, e.g.. Theorem 4.2.9 in lBleeckeii [19811): 



(^77)(d'"«,/?)/i= / {gl){a,6''P)^I, (2.4) 

AI J M 

where denotes the volume form associated to the Riemannian metric g. 

Given a connection A, a Riemannian metric g on M, and an Ad-invariant inner 
product 7 on g, we can define the Kaluza-Klein metric Kji^ on TP by 

Ka{up, Vp) := g{Tpn{up),Tpn{vp)) + -f{A{Up), A{vp)). (2.5) 

The Kaluza-Klein metric is G-invariant, that is, = Kj^ for any g E G. 
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2.5 Expressions in a local trivialization 

Consider a local trivialization Pu := 7r^^(f/) ^ U x G, p ^ {x,g). This induces a local 
trivialization of the vector bundle Ad P M, given by 

AdPDPuXGQ^UxQ, [{x,g),^]G^ {x,Adg^). (2.6) 
It is useful to note that for u G Q'^{P, g), we can locally write 

uj{x,g){{vi,ai), (ffc,afc)) = Ady-i(p{x){vi, ...,Vk)), 

where uJ{x){vi, ...,Vk) := u;(x, e)((f i, 0), {vk, 0)) e Q, Vi e TJJ, G T^G, i = 1, . . . , fc. 
Moreover, in the local trivialization fl2.6p . the Ad P- valued fc-form uj G fi'^(M, AdP) 
defined in ( 12.11) is given by 

uj{x){vi,...,Vk) = {x,uj{x){vi,...,Vk)), xeU, Vi^TJJ. 
Recall also that, in a local trivialization, a connection A can be written as 

A{x,g){v,,Q = Mg-.{A{x){v^,) + TRg-.{Q), (2.7) 

where Vx G T^IJ , G TgG^ and ^ is a 1-form on U (Z M with values in g. Locally, for 
/ G J?-g(P, g) and the curvature i3, we can write 

&^f{v)=df{v) + \A{v)Ji 
B{y,w) = dA{v,w) + [A{v),A{w)]. 

If the principal bundle is trivial, the previous formulas hold globally and the adjoint 
bundle is also trivial AdP = M x g, so we have fi''(M,AdP) = n''{M,g) and for 
u G Q'^^P, g) we get u = lj. 



2.6 Automorphisms and gauge transformations 

We say that a diffeomorphism of P is an automorphism if it is equivariant, that is, 
o ip = ip o $0, for all q E G. The Frech et Lie group of all automorphisms is denoted 



by Aut(P). See iKriegl and Michorl |1997l | for an account of Frechet Lie groups in the 
framework of manifold of maps from the point of view of the "convenient calculus" . An 
automorphism of P induces a unique diffeomorphism ^ of M defined by the condition 
n o (f = Tp o 71. The Lie algebra aut(P) consists of G-invariant vector fields on P. Its 
(left) Lie bracket is denoted by [U, V]l and is the negative of the usual Jacobi-Lie bracket 
\U, V]jL- For U G aut(P) we denote by \U] G X(M) the unique vector field on M defined 
by the condition 

T-KOU =[U]o TT. 

The subgroup Aut^{P) consists, by definition, of automorphisms ip oi P such that the 
induced diffeomorphism ip preserves the volume form /i on M. For U G out^(P) we have 
[U] G Xdiv(^), the space of all divergence free vector fields on M. 
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The normal subgroup Qau{P) of gauge transformations contains, by definition, 
all automorphisms (f on P with Tp = idM- Note that we can identify the gauge group 
Qau{P) with the group 

jr^(P, G) := {r G J^(P, G) | r o $^ = AD^-i o r}, where M)g{h) := 

The identification is given by the group isomorphism^: Qau{P) —>■ J-'g{P,G), which 
associates to ip E Qau{P), the map (p G J-'g{P, G) defined by the condition 

vip) = '^mip)- 

The Lie algebra gau(P) consists of G-invariant vertical vector fields on P. Therefore 
when U G gau(P) we have [U] = 0. Note the identifications 

0au(P)^ J-G(P,g) = r(AdP). 

Indeed, to / G J-'g{P,q) we can associate the G-invariant vertical vector field a{f) G 
Qau{P) given by 

^if)ip)--=fip)p{p)- (2.8) 
The second isomorphism is given by the map (12 .11) . A direct computation shows that 
cr : J-'g{P, q) — »■ 0Ciu(P) is a Lie algebra isomorphism, that is, 

a{[f,g]) = [aif),aig)],. 

The transformation law of a connection A under (f G Qau{P) is given by 

If* A = Ad^-i + TL^-i o T(p. (2.9) 

If the principal bundle P — > M is trivial the automorphism group is the semidirect 
product of T>{M) with J-'{M, G). To see this, note first that each ip G Aut{P) is in this 
case of the form 

where Tp G V{M), the diffeomorphism group of M, and ^ G JF(M, G), the smooth G- 
valued functions on M. Thus the map ip G Aut{P) ^^ (jp,^) G ©(M) x T{M,G) is 
bijective. Second, the pair (^r"o"^, corresponding to the product ipi o is 

uniquely determined by the right hand side of the identity 

(y^i 0(^2) {x,g) = ((^o^)(x), (^o^)(x)^(x)^) . 

This shows that the map ip G Aut{P) ^) G V{M) (§)J^{M, G) is a group isomor- 

phism, where the semidirect product is defined by the right action of V{M) by group 
automorphisms on jr(M, G) given by (x, A) G J^(M, G) x V{M) ^ x o A G J^(M, G). In 
particular, ii ip E Qau{P), then ?/' = id-M and we have ip{x,g) = (x, ■i/j{x)g), which shows 
that the map ip ^ Qau{P) ^ ip E J-'{M, G) is a group isomorphism. 

The same considerations hold for the volume preserving case. We have shown hence 
that if the principal G -bundle it : P ^ M is trivial, then we have the group isomorphisms 

Aut{P) ~ V{M) (s) gau{P) and Aut^{P) ~ V^{M) (s) gau{P) 
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and the corresponding Lie algebra isomorphisms 



aut(P) ~ X(M) (S)J^(M, g) and aut^(P) - Xaw{M)(§)T{M, g). 
Using the general formula for the Lie bracket associated the L i e alge bra of a semidirect 



product of two groups (see formula (6.4.2) in iMarsden et. al.l [20071], for example), we 



find that the (left) Lie bracket on aut(P) and aut^(P) is 

[{v, 9), (w, co)]l = i[v, w]l, d9{w) - duiv) + [9, cu]) . (2.10) 

If the principal bundle P — > M is not trivial, the situation is more involved. First, 
the sequence 

— ^ gou(P) — > aut(P) — > X(M) — > 

is exact. The second arrow is the inclusion and the third is the Lie algebra homomor- 
phism given by L'" G aui(M) [JJ] G X(M) which is surjective because any X G X(M) 
is covered by its horizontal lift relative to some connection. Note, however, that the 
horizontal lift of vector fields relative to a connection is, in general, not a Lie algebra 
homomorphism since the bracket of two horizontally lifted vector fields has a vertical 
part. This is an indication that if P — M is nontrivial, then aut(P) is, in general, not 
the semidirect product of X(M) with gou(P). 

Second, at group level, the map Aut{P) — > V{M) is not surjective, in general. For 
example, let P = S''^ C M'^ be the unit sphere, thought of as the unit quaternions, and 
let 5*^ := {x + yki x,y & M, x'^ + y'^ = 1} act on P by g i-^ q{x + yk). The Hopf fibration 
map TT : q ^ ^ qkq G S"^ defines a principal S^-bundle over M := S*^. We shall prove 
that in this case the map Aut{P) V{M) is not surjective0 Let t] G V{S'^) be the 
antipodal map whose degree is -1 and is hence not homotopic to the identity. If there 
were some ip G Aut{S'^) descending to r/, then ip would not have any fixed points and 
hence its degree would be one. By the Hopf Degree Theorem ip would then be homotopic 
to the identity which would imply that rj was homotopic to the identity, a contradiction. 



2.7 Duality 

In this paper we will identify the cotangent space T*Aut{P) with the space of G-invariant 
1-forms on P along ip G Aut{P). The duality pairing is 

(M^,[/^) := / M^(f/^)/i, 

J M 

where M,^ G T*Aut{P) and G T^Aut{P). Note that in this formula we used the fact 
that M(^([/(^) is a smooth function on P that does not depend on the fiber variables and 
hence induces a unique smooth function on M which is then integrated using the volume 
form /i on M. In particular we have aut(P)* = fiJj(P), the space of right-invariant 1- 
forms on P. 

We identify the cotangent space T*Qau{P) with the tangent space T^Qau{P) via the 
duality 

{U^,V^) := [A{U^),A{V^)) /i, (2.11) 
^We thank Marco Castrillon-Lopez for this example. 
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for any principal connection A on P. Note that, since and are vertical, the pairing 
(12. lip does not depend on A since A{Uy,) = a~^{U^ o ip~^) o for any connection A. 



3 Equations for the fields 

In this section we give the Lagrangian and Hamiltonian formulations for the Yang-Mills 
fields in the vacuum. We will see that it is not possible to pass from one to the other by 
a simple Legendre transformation. 



3.1 Lagrangian formulation of the Maxwell equations 

On the Lagrangian side, the variables are the magnetic potential A G VL^{M) and the 
electric potential Aq G JF(M), where M is a three dimensional compact manifold without 
boundary. The Lagrangian is defined on the tangent bundle T(jF(M) x fi^(M)) and is 
given by 

L{Ao,Ao,A,A) = l- [ \\Ef^-l[ \\Bffi, 

where E := —A + dAo, B := dA, and || • || is the norm associated to the Riemannian 
metric induced by g on the vector bundle K^M —>■ M, for k = 1,2. The Euler-Lagrange 
equations associated to L are 

dE 

dE = and — = 6B. (3.1) 

dt ^ ' 

The relations E = —A ^ dAg and B = dA give 

dB 

— = -dE and dB = 0. (3.2) 

dt ^ ^ 

Using the vector field variables E := and B := (*-B)^ where -k : Q''{M) — > Q^~''{M), 
is the Hodge-star operator associated to the Riemannian metric g on M, we obtain the 
Maxwell equations in the vacuum 

dE 

— = curlB, divE = 0, 

— = -curlE, divB = 0, 
at 

where curl : X(M) X(M) is the operator curl(X) := [7^(dX^)]f, for any X G X(M). 

Let us recall the classical argument that we can choose Aq = 0. Assume that Aq 
and A' satisfy Maxwell's equations ( 13. ip and ( 13. 2p . We search a function ip G JF(M) 
such that A := A' + dip leaves the equations (13. ip and (13. 2p unchanged and Aq = 0. 
Since these equations are second order we have A'o = dAQ/dt and A' = OA' /dt. Let 
E' := —A' + dA'Q, B' := dA'. The requirement is that E' = E and B' = B. Therefore, 

-A' + dA'o = E' = E = -A + dAo = -A' - d0 + d^o 
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which is equivalent to dip = d(ylo — Aq) and hence it is sufficient to choose Aq := A^ + ip. 
This shows that one can choose Aq = provided (p = —Aq. Note that the equations 
(13. ip and ( 13. 2p are unchanged under this transformation, as required. 

We now recall the four dimensional formulation of the Maxwell equations. Consider 
the Lorentzian manifold (X, 7) given by X = M x M and 7 := r^g — r2(it^, where 
Ti : X ^ M and T2 : X —>■ M. are the natural projections and g is a Riemannian metric 
on M. Consider the 1-form G on X defined by G := T*At + t*Ao A rlc/t. We have 

dG = T*dAt - T*At A T^dt + dAo A r^dt = t*B + t*E A r^dt =: F, 



and th e Maxwell equations (13. 3p can be simply written as (see e.g. iMisner. Thorne. and Wheeler 
1973j . §22.4) 

dF = and 5F = 0. 
In a general slicing of space-time, not just M x M, the de rivation of these equations and 



much more information can be found, for example, in iGotay. Isenberg. and Marsden 



1997, 1999 



The Legendre transformation associated to the Maxwell Lagrangian L is not bijec- 
tive. Thus, it is not possible to pass in the usual way from the Lagrangian to the 
Hamiltonian formulation by the Legendre transformation. This degeneracy is typical 
of relativ istic field theories and is resolved by the D irac theory of constraints; see, for 
example, iGotay. Isenberg. and MarsdenI 19971 . Il999l | and references therein. In the next 
subsection we directly generate the Hamiltonian formulation for the Maxwell equations. 



3.2 Hamiltonian formulation of the Maxwell equations 



On the Hamiltonian side (see lMarsden et. al.l 1983| ). the configuration space variable is 
the magnetic potential A e Q^{M). The Hamiltonian is defined on the cotangent bundle 
T*Q^{M) ~ Q^{M) X Q^{M), where the cotangent space at any point A is identified 
with Q^{M) using the natural L^-pairing, and is given by 



J M J M 



6B, 



for E := —Y and B := dA. Hamilton's equations are 
and the relation B = dA gives 

dB = 0. 

To obtain the last equation 6E = we use the invariance of the Hamiltonian under 
gauge transformations. The action of the gauge group JF(M) on Q^{M) is given by 



J^{M) xn\M) ^Q\M), {^,A)^A + d<p, 
and is Hamiltonian. The associated momentum map is 

J : T*fi^(M) ^ J^(M)* ~ J^(M), J{A,Y) = 5Y, 



(3.4) 
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where JF(M)* is identified with J^{M) using the natural L^-pairing. So the condition 
J {A, Y) = gives the fourth Maxwell equation 6E = 0. 

Note that in the Hamiltonian formulation we have used only the configuration vari- 
able A, whereas in the Lagrangian formulation the configuration space consisted of pairs 
{Ao, A). As we have seen, the variable Ao can be set equal to zero without any effect on 
Maxwell's equations. Note also that the Euler-Lagrange equation 6E = was obtained 
from the variation of the Lagrangian relative to Aq, whereas in the Hamiltonian set-up 
this equation appears as a conservation law for the gauge group action (13 ■4p . 

3.3 Generalization to any principal bundle 

We now generalize the previous formulations to the case of a G-principal bundle P M 
over an arbitrary compact boundaryless manifold M. We will show that if M is three 
dimensional, G = S^, and the bundle is trivial, then we recover the Maxwell equations. 

Lagrangian formulation. The Lagrangian L : T{J-'g{P, q) x Conn[P)) ^ M is defined 

by 

L{Ao,A,,AA) = \ f WEf^^-l [ (3.5) 

^ Jm ^ J M 

where: 

(1) E := S & Q^{M,AdP) is the AdP-valued 1-form associated, through the map 
f lO) . to the "electric part" £ G n^{P,g) of the Yang-Mills field, given by 

S := -A + d^AoeW{P,gy, 

(2) B := B E f2^(M, AdP) is the AdP-valued 2-form associated, through the map 
(12.11) . to the "magnetic part" B E f2^(P, g) of the Yang- Mills field, given by the 
curvature 

B:=d^Ae'D?{P,Qy, 

(3) II ■ II is the norm associated to the metric (77 on the vector bundles A^{M, Ad P) — ^ 
M, for k = l,2. 

The Euler-Lagrange equations associated to L are 

df 

S^£ = and — + [Aq, £] = 5^B. 
ot 

Indeed, using the pairing 

(«,/?)=/ {gi){a,P)fi, a,Pen\M,AdP), (3.6) 

we can identify the cotangent bundles of J^g{P,q) and Conn{P) with their tangent 
bundles. Using formulas (12. 3p . (12. 4p . and the identity 

d^+tc^A + tC) =d^C, (3.7) 



d 
di 
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where A G Conn{P), C G TjjZonn^P) = fi^(P, g), we get 

Thus, the Euler-Lagrange equations 

d dL dL ^ d dL dL 

t; ■■ TT-^ = and — — — = 

dtdAo dAo dtdA dA 

become 

d£ 

5-^8 = and — + Uo, ^] = 5"^^, 
dt 

as stated above. 

The relations £ := — ^ + d"^^o and B := d"^^ give the equations 
dB 

+ Uo, B] = -d^S and d-^B = 0. 

at 

Indeed, for the first equahty we have, using (13.71) and Lemma [2TT| 

B = d-^A = -d-^S + d-^d-^^o = + [B, A]- 

The second equahty is the Bianchi identity. Summarizing, we get the system 

df 

— + [A, £] = 5^B, = 0, 
dB 

— + [Ao,B] = -d^£, d^B = Q. 



(3.8) 



To recover Maxwell's equations we take a trivial 5'^-principal bundle P = M x S^. 
Then AdP = M x R and n''{M,AdP) = Vt^{M). Since the structure group of 
the principal bundle P is Abelian, the covariant differential does not depend on the 
connection, that is, d"^ = d. We obtain the following identifications. 

(1) From the equality £ = —A + d'^Ao, we obtain that the electric field E := £ G 
il^(M) is given by 

E = -A + dA, 

where A := A e J^(M) and A := A e Q^{M). 

(2) From the equality B = d"^^, we obtain that the magnetic field B := B E f2^(M) 
is given by 

B = dA, 

where A G Q^{M) is given hj A := A (see equation (12.71) ). 
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Returning to the general case, let us show, as in the case of Maxwell's equations, 
that we can choose = 0. Assume that and A' satisfy equations (13.81) . We search 
a e Qau{P) such that A := (p*A' leaves the equations (13. Sp unchanged and = 0. 
Since these equations are second order we have Aq = dA'Q/dt and A' = dA' /dt. Let 
£' ■= -A' + d-^'A'o, B' := d^A'. Since B = ip*B', the requirement is that £ = ip*£'. 
Therefore, emphasizing the time- dependence, we have the equivalences 



At + d^'Aot = ^;[-A + d^'^A', 
di 



ot 



(3.9) 



-i^:A[) + d-^^Aot 



^,-A[ + d^^^;X, 



t=s 



Taking the time derivative of (12. 9p we get 

rtA = d^{a-\^o: 



d 

dt 



t=0 



(3.10) 



for any smooth curve ipt £ Qau{P) such that %Ijq = id. Therefore we conclude that 



d_ 

ds 



* /I ' 



d_ 

ds 



s=t 



So (13.91) is equivalent to d'^'^ot = d"^Vt ""^(^t ° Vt'^) + At) hence it is sufficient 
to choose Aot '■= {'^~^i.'Pt° Vt^) +-^0*) order to get £ = {p*£'. Thus one can 
choose Aot = provided that 0^ o {p~^ = — o"(^qJ. A direct computation shows that the 
equations (13. 8p are unchanged under this transformation, as required. 

Let (X = M X R, 7) be the Lorentzian manifold given in §3.11 Let P := P x M. 
and define the free G-action $p(p, t) := {^g{p),t) on P. We get the principal G-bundle 
W : P := P X M. ^ X . From At € Conn{P) and ^0 ^ ^g{PiQ)i we can construct the 
1-form C e fiHP,0) 

C := T*At + T*Ao A {t2 o wydt, 

where r : P ^ P is the natural projection and r2 : X ^ M is the projection on the 
second factor. One can check that C is a connection on P since t*A G Conn{P) and 
T*Ao A {t2 owydt eW(P, q). 
Finally we obtain 

d^C = T*B + T*£ A {t2 o wydt =: 

and equations fl3.8p are equivalent t o the Yang-Mills equations together wi th the Bianchi 
identity (see e.g. lArmd 19791 . Il98l| . lArms. Marsden. and Moncriefl 1982| ) 



b^T = and d^J^ = 0. 
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Hamiltonian formulation. As in tlie electromagnetic case, the configuration space 
variable is the magnetic potential A G Conn{P) and the Hamiltonian is defined on the 
cotangent bundle T*Conn{P) by 

HiA,y) = l f WEfi^+l [ ||Sf/i, 

^ J M J M 

where: 

(1) E := S E Q^{M,AdP) is the AdP-valued 1-form associated, through the map 
f lO) . to the "electric part" £ G '^{P,g) of the Yang- Mills field, given by 

^:=-3^GlF(P,g), 



(2) B := B E f2^(M, AdP) is the AdP-valued 2-form associated, through the map 
(12.11) . to the "magnetic part" B G ^^(P, g) of the Yang- Mills field, given by the 
curvature 

B:=d^AE'Q^iP,Q). 

As before, we identify the cotangent bundle of Conn{P) with the tangent bundle, 
using the pairing fl3.6l) . 

Hamilton's equations associated to H are 

— = -d-^S and — = 6-^B, 
dt dt 

and the Bianchi identity gives 

d-^B = 0. 

To obtain the last equation we use the invariance of the Hamiltonian under gauge trans- 
formations. The action of G Qau{P) on ^ G Conn{P) is if* A and the cotangent lift 
of this action is {ip*A,(p*y). Under this action, S and B are transformed into ip*S and 
(p*B, so H is gauge-invariant. The momentum mapping associated to this Hamiltonian 
action is 

J : T*Conn{P) gau(P)* ~ gau(P), J(A 3^) = a{5^y), 
so the conservation law i{A,y) = gives the last equation 

6^S = 0. 



Note that we identify gau(P)* with gau(P) via the pairing (13. 6p . 



4 Equations for the particles 

We consider the evolution of a non relativistic Yang-Mills charged particle of mass m in 
a given Yang-Mills field. 

Fix a connection A G Conn{P) and an equivariant function ^ ^G{P,d)- The 
Yang-Mills field is given by its electric part S = d'^^o and its magnetic part B = d'^A. 
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Consider the right-invariant Lagrangian L : TP — > M, given by 



Tfi 1 

L{up) = —g[p]{T7T{up),T7T{up)) + --f{A{up) + Ao{p), A{up) + Ao{p)). 

Since L is G-invariant, it induces a Lagrangian on (TP)/G. We use the identification of 
(TP ) I G witli TM ® Ad P tfirough tfie c onnection dependent vector bundle isomorpliism 
(see 



Cendra. Marsden. and Ratiu 200lj |) 



: {TP)/G ^TM® AdP, ^^(M) := {Tpn{up), [p,A{up)]g). 
Tlie reduced Lagrangian on TM © Ad P is given by 



TTl 1 



wfie re An G FfAdP) is associated to via tlie map (12.11) . By Lagrangian reduction 



see 



Cendra. Marsden. and Ratiu |200l| ). p{t) G P is a solution of the Euler-Lagrange 
equations for L if and only if x{t) := Tr{p(t)) G M and C,(t) := [p(t) , A{p(t))]G G AdPj;(t) 
are solutions of the Lagrange-Poincare equations 

fa/,.,, D9 dl / dl 

D-^ dl dl 

where D^/dt and D-^/dt denote the covariant derivatives induced by g on T*M and by 
A on (AdP)*, respectively, 

dl dl 

g^ivx,Q e T;M and ^{v,,Q G (AdP): 

are usual fiber derivatives of / in the vector bundles TM and Ad P, and 

dl 

-{v^,QeT:M 

is the partial covariant derivative of / r elative to the Levi-Civita connection on M and 
the principal connection A on P. See ICendra. Marsden. and Ratiul 2001| for details 
regarding the Lagrange-Poincare equations. 
In terms of the functional derivatives 



J^iVx,^x), ^ivx,Q^T^M, 



and 



61 



,v,,^,) G (AdP), 



defined similarly, the Lagrange-Poincare equations become 



r 61 



D3 61 



6x 
D-^ 61 



dt 6v 



61 
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where [■, ■]x is the bracket of elements in (AdP)^ and and ^ denote the covariant 
derivatives on TM and Ad P, respectively. Using that 



where 



61 

TT(fx,^x) = + Ao{x) 



we obtain that the trajectory x{t) := n{p(t)) G M and the charge 

q{t) = m,Apit))+MpmG = ^{i{t),m) e (AdP),.(,), 



are solutions of 

^^iW = 7x(t) E{-) + B{; x{t))Y 

— g(t) + [Ao(x(t)),g(t)],(t) = 0. 

The first line is the nonabelian Lorentz equation and the second line represents the co- 
vari ant conservat i on of charge. These equations are the non-relativistic Wo ng equations 
(see 



Wond |l970l |. lMontgomervl |l984l | • ICendra. Marsden. and Ratiul |200l| l 



In the case of the trivial S'^-bundle P = M x S'^, recall that AdP = M x M. The 
Lagrangian is 



Tfl 1 

-^xK, v^) + -{A{x){v,) + 6 + Ao{x)f, 



where A := A & n\M) and Aq := Aq. We have {TP)/S^ = TM ®MP = TM x R, 
and 0) = {vx, A{x){vx) + 0). So the reduced Lagrangian is 

Tfl 1 

^K, = —9x{Vx, Vx) + -(^ + Ao{x)Y. 



By Lagrange- Poincare reduction f Cendra. Marsden. and Ratiu 200lj |). we obtain that 
pit) = {x(t), 9(t)) is a solution of the Euler-Lagrange equations if and only if 

m!^x(t)=q(t)(E(-) + B(;x(t)))l 
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where q{t) := A{x(t)) + 9(t) + Ao(x(t)) is the charge. If diniM = 3, in terms of the 
vector fields E := and B = (^-B)" and using that the charge q(t) = q is conserved, 
the previous system becomes simply the Lorentz force law 

m—x{t) = q{E + x{t) x B), 

describing the motion of a charged particle of mass m in a fixed electromagnetic field. 

We remark that when the variable is absent, the Lagrangian is given by the 
Kaluza-Klein metric, 

1 TTl 

L{up) = -Kj^{p){up,Up) = —gip]{Tp7i{up),Tnp{up))+-f{A{up),A{up)). 

In this case the Legendre transformation is invertible and the corresponding Hamiltonian 
on T*P is ^ 

H{ap) = -K;^{p){ap,ap), 
where K*^ is the dual metric on T*P, defined by 

^liP) iKA{p){up, ■),Kj^{p){vp, ■)) := Kj^{p){up,Vp). 



5 Lagrangian formulation of Euler- Yang-Mills 

We begi n by recalling some facts about E u ler-Poincare reduction for semi direc t prod - 
ucts (see 



Holm. Marsden. and Ratiul |1998l |. iCendra. Holm. Hoyle. and MarsdenI [1998 



Cendra. Marsden. and Ratiul 2001| ). Let p : G ^ Aut(K) denote a right Lie group 
representation of G in the vector space V. As a set, the semidirect product S = G(§)V 
is the Cartesian product S = G x V whose group multiplication is given by 

{9uVi){g2,V2) = {gig2,V2 + pg^ivi)). 

The Lie algebra of 5* is the semidirect product Lie-algebra, s = 0©^) whose bracket 
has the expression 

ad(5,,^,)(6,f2) = [(6, ^^i), (6,^2)] = ([6, 6], ^16 -^2-^1), 
where denotes the induced action of g on V, that is. 



, d 



t=0 



From the expression for the Lie bracket, it follows that for v) E s and (/x, a) G s* we 
have 

ad*g a) = (ad^ p + voa,a^), 
where G V* and v o a E Q* are given, respectively, by 



where (-,■)„ : Q* x g 



dt 



t=Q 

and {■, ■)y : V* X V ^ M. are the duality parings. 
Lagrangian semidirect product theory with parameter. 
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• Let Q be a manifold on which G acts trivially and assume that we have a function 
L :TG xTQ xV* which is right G-invariant. 

• In particular, if & V*, define the Lagrangian L^;, : TGxTQ — M by La„{t-'g: Uq) : = 
L{vg,Uq,aQ). Then La^ is right invariant under the lift to TG x TQ of the right 
action of Gao on G x Q, where Gao is the isotropy group of gq. 

• Right G-invariance of L permits us to define I : QX TQ x V* ^ R hy 

KTgRg-^{Vg),Uq,p*g{ao)) = L{Vg,Ug,ao). 

• For a curve g{t) G G, let C,(t) := TRg(^t)-^{9{t)) and define the curve a(t) as the 
unique solution of the linear differential equation with time dependent coefficients 
a(t) = —a(t)^(t) with initial condition a(0) = oq. Its solution can be written as 

= P*g{t)i(^o)- 

Theorem 5.1 The following are equivalent: 
i Hamilton's variational principle holds: 

d r LaMt).g{t),q{t)A{t))dt^Q, 

for variations of g and q with fixed endpoints. 

ii {g{t),q{t)) satisfies the Euler- Lagrange equations for La^ on G x Q. 

iii The constrained variational principle 

S mt),Qit),q{t),a{t))dt = 0, 
Jti 

holds on QX Q, upon using variations of the form 

Oti 

0^ -[^^V], Sa^-ar], 

where rj{t) G Q vanishes at the endpoints and Sq{t) is unrestricted except for van- 
ishing at the endpoints. 

iv The following system of Euler- Poincare equations (with a parameter) coupled with 
Euler- Lagrange equations holds on Q x TQ x V*: 

d 51 ^,51 51 

and 

d_dl__dl_^ 
dt dq dq 
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Note that the Euler-Poincare equation (15. ip can be written, in weak form, as 



-D((0(K,.;l) + (^oa.,, 



for all rj E g, 



(5.2) 



where D denotes the Frechet derivative. This formulation will be useful below. 

Ideal compressible adiabatic fluid. Before treating the Yang-Mills fluid, we apply 
the preceding theory to the case of the compressible adiabatic fluid. For this particular 
case we choose G = V{M) and V = JF(M) x JF(M) (in this case Q is absent). We 
identify the dual JF(M)* with J-'{M) via the natural pairing. The action ofrjE T>{M) 
on (p, s) e V* is 

(p, s) ^ {{Jr]){por]),sori), 

where Jrj is the Jacobian determinant of rj, p is the density of the fluid and s is its specific 
entropy. As usual, we treat the mass density p as a density on M and the entropy s as 
a function on M; this is why in the previous formula the action of the diffeomorphism 
group is different on the two components. 
The Lagrangian is given by 



1 



■'{p,s) \ "'Tj 



M 



pg{urj,Ur,)fi- / pe{p{Jri) \ s)p. 



(5.3) 



M 



where e is the fluid's specific internal energy. Application of part iv in Theorem 15.11 
gives the equations of motion 



dv „ 
-7^ + V^,f 



dt 
dp 
dt 



— gradp, 
P 



+ div(pf ) = 0, 



(5.4) 



ds , , , 
de 



where the pressure is given by p = p — (p, s). 

op 



5.1 Yang-Mills ideal fluid 

In the case of the Yang-Mills fluid we choose G = Aut{P), Q = J^ciP^d) ^ Conn{P) 
and V = J^{M) x J^{M). As before, we use the notations (f G Aut{P), {Ao,A) G 
J-'ciP, q) X Conn{P), and (p, s) & V*. The action of ip on (p, s) is given by 

(p,s) ^ ((J^)(po^),so^), 

where ^ G T>{M) is the map induced on the base M by ^p. From the expressions of the 
Lagrangian (15.31) and of the Lagrangians for the fields and particles given in ^ and ^ 
it follows that the Lagrangian for the Yang-Mills ideal fluid is defined on the tangent 
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bundle T{Aut{P) x ^"^(^,0) x Conn{P)) by 



L^p,,){U^,Ao,Ao,A,A) = 1 I pgi[U^],[U^])fi (5.5) 

^ 'M 



+ \ p\\{A{U^)+Aqo^) II V 
- / pe(p(J?/')"\s)/i 

J M 

+ 1 /■ II^IIV-- / ||5f/i, 

where [U^] e T-:^V{M) is such that Tn o = [U^] o tt. Note that A{U^) + Aq o e 
J-'g{P, 0), so we can consider the section {A{U^) + ^0 ° ^)~£ r(AdP) and its norm 
II {A{U^) + ^0 ° relative to the Riemannian metric 'jx- The two last terms of (15.50 
are given as in the Lagrangian (13. 5p . Roughly speaking, this Lagrangian has the following 
structure 



Integration of the Lagrangian 
for the particles 

Lagrangian for the 
Yang-Mills fields 



{Internal energy} 



+ 



Note that L verifies the invariance property needed for an application of Theorem 
15.11 that is, L is invariant under the right action of G Aut{P) 

{U^,p,s) ^ {U^oip,{Jlp){polp),soJp). 

Indeed, we have [U^oip] = [t/^] o^, so the invariance of the first term follows by a change 
of variable in the integral. The invariance of the second integral follows from the fact 
that A{U^) + Aqo^P e Tg{P-, 0) and that for /, (7 G Tg{P-, 0) and ^9 G Aut{P) we have 

7(/ o 9°^) = 7(7, 9) o ^, 

as functions on M. 

The reduced Lagrangian / on aut(P) x T{J^g{P^ 0) x Conn{P)) x (J^(M)* x J^(M)*) 
has the expression 

l{U,Ao,Ao,AAp,s) = \ f pg{[UUU])p (5.6) 

+ \ pUA{U)+Aorfp - I pe{p,s)p 

^ JM J M 

+ 11 \\E\?p- \ I m^p 

^ JM ^ JM 

and the Euler-Poincare equations in weak form are 

|D/(f/)(F) = -D/(f/)([[/, \/]^) + o (p, .), , (5.7) 
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for all V G aut(P). We now compute these equations. 

Recall that the (left) Lie bracket on the Lie algebra out(P) is 

[U,V]L^SiduV = -[U,VU, 

where [ , ]jl denotes the usual Jacobi-Lie bracket of vector fields. The following lemma 
gives the decomposition of [U, V]l into the horizontal and vertical parts. 

Lemma 5.2 Let A he a connection on the principal bundle P and let U,V & out(P). 
Then we have 

[U,V]l = a{[A{U),A{V)]+d-^iA{U)){V) (5.8) 
- d^{A{V)){U) + B{U, V)) + hor([t/, VU), 

where hor denotes the horizontal part relative to the connection A. In particular we have 
the equality 

A i[U, VU) = [AiU),AiV)] + d-^{AiUmV) - d-^(^(V))(^) + BiU, V). 

Proof. First note that using the Cartan Structure Equations and the fact that B e 
02(P, 0), we have 

dA{hor U,a{A{V))) = B{horU,a{A{V))) - [A{horU),A{a{A{V)))] = 0. 
We also have 

d^(horC/,(7(^(y))) = d(^((7(^(y))))(horC/) - d(^(horC/))((7(^(y))) 

-A{[horU,a{A{V))]jL) 
= d(^(y))(hor U) + A{[hoi U, a{A{V))]L). 

These formulas prove that 

d{A{V)){horU) = -A{[horU,a{A{V))]L). (5.9) 

We now compute the Lie bracket [U, V]^. By decomposing U and V into their vertical 
and horizontal parts, that is, we write U — a{A{U)) + hor U and V — a{A{V)) + hor V, 
we obtain four terms. The first term is 

[hor [/, hor V]l = cr (.4. ([hor U, hor V]l)) + hor [hor U, hor V]l 
^a{B{U,V))+hor[U,V]L, 

where we used the equalities 

B{U, V) = -^([hor U, hor V]jl) 

and 

T7ro[U,V]L = m,[V]]L0 7r. 
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Since [a {A{U)), hoi V]l is vertical (apply the formula above), the second term is 

[a{A{U)),horV]L = aiAMAiU)),hoTV]L)) = d(^(f/))(hor V) 
= d-'{A{U)){V), 

by formula (15.91) . There is an analogous formula for the third term [hoTU,a{A{V))]L. 
Using the Lie algebra isomorphism a : jFc.(P, g) Qau{P) defined in (12. 8p . the fourth 
term is 

[a{A{U)),a{A{V))]=a{[A{U),A{V)U). 

Summing these four terms we obtain the desired formula (15. 8p . ■ 

Inspired by the Kaluza-Klein metric (12.51) . we define on aut(P) a non-degenerate 
bilinear form given by 



Therefore we have 



M 



g{[U],[V])fi+ [ j(AiU),AiV))fi. 



^W,\U,V\l)j, = [ g{[W],[[U,V]L])fi+ [ ^(A(W),A{p^^]L))l.i (5.10) 

g{\Wl [[t/], [\/]]^)/z + / 7 [m\ \A{\J\A{V)^ 
M Jm ^ ^ 



- 7 



M \ J 



+ / 7(^),S([f/],[^]))/i 
Jm ^ ' 

+ 7 (-Aof), A\Mum + mi ■)) \\y\]i^ 
+ / ^{\A{v^m-.d^Xm)w\ 

Jm \ 

+ div([t/]).A(HO, ^1 



where in the last equality, ad^ denotes the adjoint of ad relative to the metric and 
jj is the index raising operator associated to g. Note that for m, w e X(M), ad'^ IS given 
by 

ad|^ w = V„w + Vu"^ ■ + If div u. (5-11) 
In the second summand of the last equality in ( I5.10p we used the following Lemma. 



24 



Lemma 5.3 Consider an Ad-invariant inner product 7 on g and the induced vector 
bundle metric on AdP, also denoted by 7. Then for v G X(M) and f,g& J-'g{P,q) we 
have 

(d^f{v),g) fx = - (f,d^g{v)) fx- j{f,g){divv)fi. 

^ ^ JM ^ ' JM 



' M 

Proof. One verifies that for any v G X(M) we have 



d {1U.9)) {v) = 7 {d^f{v),g) +^[f,d^g{v: 



Integrating over M gives the result. Indeed, denoting by h the real valued function 
7(/, g), we obtain 



h{divv)fi, 



/ dh{v)fi= / div(/if)yU— / h{divv)ij, = — 
J M J M J M 

by the divergence Theorem. ■ 

Using flS.lOp and the formula 

we obtain 

D/(f/)([f/,y]i) = (pf/ + aMo),[f/,^]L)^ (5.12) 
= g l^adf^j p\U\ + 7 (q, d-V(f/))(-) + ^([f/], •)) " , \y\\j 

[qT^)] + d^Q{V) + div([[/])Q,J(V)) /X, 



+ / 7 

J M 



where 



Q := p{A{U) + A) G g) (5.13) 

is the charge density. 

On the other hand we have, using the notations Aq := Aq E r(AdP), A := A E 
n^{M,AdP), and Q := Q E r(AdP), 

= 9 {^/P], [V]^ p + j^p (q, MY)) p + j^ i (q, A[V]) p 

= j^^g (p|[f/] - div(p[[/])[t/] + 7 {Q.K))\ [V]^ P 

+ / 7(43oo)p, 

J M ^ ' 
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where we used the equation p = — div(p[t/]). Using the equahties 



Q Q 
P ' P 



e- p 



de 
dp'' 



6l_ 
Ts 



-P 



de 
ds' 



51 , , 51 61 ^ 

^^-^o(p, s) = pgrad — - — grads, 

1 



pV[UY ■ [U] + -pgrad7 (-,-)- grad ( p' — 



Q Q 



P P 



,de 



equation (15. 7p yields the system 



d 



p^[f/] + pV[^] [f/] + 7 ( Q, A{.) + d^{A{Um + Bi[U], 



-p grad 7 



[Q Q\ _ 
\P' p) 



gradp 



+ [Q, AiU)] + d^QiU) + div([f/])Q = 0. 
Q 



Denoting S := — = A(U) + several applications of Lemma 1531 give 
P 



M 



-f[S,d-^{A{U)){v)]p 



M 



M 



d{^{S,S)){v)p- I ^[S,d^Ao{v 



for all V G X(M). So we get 

7 (^g,d-V(t/))(-)y = ^pgrad7 ^) - 7 

With this formula and the equality 

[Q,AiU)] = [Ao,Ql 

the system (I5.14p is equivalent to 

^ [U] + V[^][f/] = -7 (q, -i(-) + f?A(-) + B{; [U]) 



dt 



d 



P 



Q + [A, Q] + d^QiU) + div([f/])Q = 0, 



— gradp 
P 



\ dt 

which is the same as 



-^ + V,v = -7 (g, E{-) + v)f - - gradp 
ot p p 

^ + [Ao,g] + v;fg + gdivt; = o, 
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where v := [U] & X(M) is the Eulerian velocity. 

We compute now the Euler-Lagrange equations relative to the Lagrangian (15.61) and 
the variables [Aq, A). The computations are similar to those done in Paragraph 3.3. We 
find 

where Q (g) Tr*v^ G (P, q) is given by 

(Q ® 7r*t;^) (up) := (77*^;^^)) = Q{p)gM^),Tp7i{up)), x = 7i{p). 

Note that we have Q vr*f'' = Q ® f'', the 1-form on M, with values in AdP, given by 

[Q^v^] (m,) = Q{x)v\u,) = Q{x)g,{v{x),u,) e (AdP),. 
For the computation of the partial derivative we use the identity 

^{Q,A{v)) = {g^){Q®v\A). 

The resulting Euler-Lagrange equations are 

BE 

5^8 = -Q and — + [A, S] = 5^B-Q® ti*v\ 
ot 



As in Paragraph 3.3, the relations £ = —A + d"^^o s-nd B = d'^A give the equations 

SB 

— + Uo, B] = -d^S and d^B = 0. 
ot 

Summarizing, we have proved the following theorem, which is one of the main results 
of this paper. 

Theorem 5.4 Let {iIj,Ao,A) be a curve in Aut{P) x jFfj(P, g) x Conn{P) and consider 
the curve {U, Ao, A) := [tpoip-'^ , Aq, A) m aut(P) x:Fg'(P, g) xConn(P). Then [ip, Aq, A) 
is a solution of the Euler-Lagrange equations associated to the Lagrangian ^v(p„^s,)) given 
in ( 15. Sp if and only if {U, Aq, A) is a solution of the Euler- Yang-Mills equations: 

( dv 1 « 1 

_ + V.i; = -7 (Q, + P(-, vyf - - gradp, 

do ds 

+ div(/w) = 0, p(0) = po, ^ + ds(^;) = 0, s(0) = Sq, 



^ + [Ao, g] + v;:^g + Q div 1; = 0, 

d£_ 
'dt 
dB 
'dt 



(5.17) 



+ [Ao,£] = S-^B-Q®7r*v\ 5-^8 



+ [A, = -d-^^, d-^B = 
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where 



P 

S 
B 
Q 



; := [U] e X(M), 
-A + d-^Ao eW{P,g) and E :=£ eVl\M,MP), 



d-^A and B := B, 
piAiU) + Ao) and Q = Q. 

Corollary 5.5 In the case of the trivial bundle P = M x and assuming that the fluid 
is composed of particles of mass m and charge q, we obtain the Euler-Maxwell equations 

dv 



q 1 

+ V-uV = — (E + wxB) gradp, 

m p 



dt 

^ + div(/w) = 0, p(0) 



(9E _ 
'dt ~ 
divE 



curl B 
(1 



m 



q 

—pv, 
m 

divB 



= Po, 

OB 

'dt 

0, 



^ + ds{v)=0,s{0) 
- curl E, 



(5.18) 



where 



E := and B := 



Proof. If we define Qt = pt^, the equation for Q in f l5.17p becomes 

= ^ + ^^t{v) = j^qt{x{t)), 

where x{t) is the trajectory of the particle starting at a;(0). Since all particles have the 
same charge g G M by hypothesis, we conclude that q(t,x) is a constant. Therefore, the 
equation for Q in fl5.17p disappears. It is easily seen that the other equations become 
the ones in (15.181) . ■ 

We end this section by examining more carefully the case of a trivial principal bundle 
P = M X G. We use the fact, already pointed out in the introduction, that in this case 
the automorphism group is a semidirect product of two groups. 

In the trivial bundle case, we have a connection independent pairing on aut(P), 
given by 



g{m,v)fi+ / -f{u,9)n 



At 



M 



Using this pairing, the expression (I2.10p for the Lie bracket on the semidirect product 
Lie algebra, the expression (15. lip , and integration by parts, we obtain the following 
expression for ad^: 



ad{^ g)(m, z/) 



(5.19) 



(ad];m + 7(z/, d6'(-))^, z/divi; + dz/(t;) + [u, 9]) 
(V„m + Vf'^ ■ m + m divv + 7(1^, d6'(-))^, 
u div V + diy(v) + [u, 6]) . 
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The reduced Lagrangian 

I : aut(P) X T(J^g(P, 0) X Conn{P)) x (:r(M)* x :r(M)*) 



IS 



l{v,e,AQ,Ao,AXp,s) = \ I pg{v,v)fi + l [ p\\A{v) + 9 + AoW^p 

-I pe{p,s)f, + l [ II^IIV (5.20) 

and we have 

^ = p^v + ^(A{v) + e + Ao,A{-)f) and ^^=p(A{v) + e + A,). 

The Euler-Poincare equations are 

d (51 5l\_ X (51 5l\ 51 

m[5^'5e)~~ ^"^M) [s^'sej^ ^ 

and a long direct computation gives, as expected, the system 

^ + V,^; = -7 (g, E{-) + ^;))« - - gradp, 
01 p p 

do — 

+ [^(t;) + A, Q] + dQ(T;) + Q div ^; = 0, 



(5.21) 



where 



A/ 

Q:= — =p {A{v) + e + Ao)e T{M, g), E e ^^(M, g), 5 e fi^^M, 0). 
5.2 The incompressible and homogeneous case 

In the incompressible case we choose G = Aut^{P), the Lie group of all automorphisms 
ip e Aut{P) such that Ip e T>i^[M). Since the fluid is homogeneous, the advected 
variables p and s are absent. Therefore, we can use standard Euler-Poincare reduction 
with parameters (^0,^) £ ^g{P,q) x Conn{P) (take y = in the semidirect theory). 
The Lagrangian for the incompressible homogeneous Yang-Mills ideal fluid is deflned on 
the tangent bundle T{Aut^{P) x Tg{P,q) x Conn{P)) and is given by 

L{U^,Ao,AoAA) = ]- I g{[U^UU^])p (5.22) 

^ Jm 

+ / \\{A{U^)+Aootljyfp 

^ JM 

+ 1 / \\Efp-\ [ \\B\W 
^ Jm ^ JM 
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The computations of the Euler-Poincare equations are similar to those done in the 
compressible case, except that we have div([?7]) = and we must replace formula (15.120 
by formula 

BI{U){[U,V]l) 

g I Pe ( adf^j[t/] + 7 (q, d-XiCum.) + mU], .)) ],[V]\^ 



M 
+ 



M 



7 [Q,^(t/)]+d-^Q(f/),^(y) 



where Pe : X(M) Xdiv(Af) is the projector associated to the orthogonal Hodge 
decomposition 

X(M) = Xdiv(M) © grad(J^(M)). 
We finally get the following result. 

Theorem 5.6 Let {%Ij,Ao,A) he a curve in Aut^{P) x JF(^(P, g) x Conn{P) and con- 
sider the curve {U,Ao,A) := {i> o ijj~'^,Ao,A) in aut^(P) x {Fg{P,q) x Conn{P). Then 
{ijj,Ao,A) is a solution of the Euler- Lagrange equations associated to the Lagrangian 
f l5.22p if and only if {U,Ao,A) is a solution of the incompressible homogeneous Euler- 
Yang-Mills equations: 

^ + V.t; = 7 {Q, E{-) + B{-, v)f - gradp, 



^ + [Ao, Q] + V^Q = 0, 
dE 

— + [Ao, S] = 6^B-Q® iv*v\ 6^8 = -Q, 
ot 

913 , , 



(5.23) 



d^£, d-^B = 



A, 



where 



V 

S 
B 
Q 



[U] e Xdiv(M), 

-A + d-^Ao eIl^{P,g) and E :=£ eVt\M,MP), 
d-^A and B := B, 
A{U) + Ao and Q = Q. 

Note that the pressure is in this case determined from v, E, and B through the 
Green's function of the Laplacian on M. This is in contrast to (I5.17P where the pressure 
was given by the internal energy. 

If P is a trivial bundle, one gets the incompressible homogeneous version of the 
Euler- Yang-Mills equations (corresponding to the group P^(M) (s) JF(M, G)) by replac- 
ing in formula ( 15.190 the vector fields by their projection onto their divergence free part, 
namely, 

ad^(m, u) = (adt m + P, d6{-)f) , du{v) + [u, 9]) 

= (Pe (V„m + Vw^ ■ m + 7(z/, d^(-))") , diy{v) + [u, 6]) . (5.24) 
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One can also adapt our method to the case of the incompressible but non-homogeneous 
Yang-Mills fluid. It suffices to apply the semidirect product theory with G = Aut^{P), 
Q = J-'g{P, 0) X Conn{P) and V = J-'{M), where the mass density p is an element of V*. 
Note that in geophysical incompressible fluid dynami cs there is also a second scalar ad- 



vect ed quantity, namely the buoyancy (for details, see iHolm. Marsden. and Ratid |l998 



and Holm. Marsden. and Ratiu 2002| ) which plays the role that entropy plays in a com- 
pressible fluid. In this case we would take V = T{M) x JF(M), where the second factor 
is thought of as the space of densities on M, thereby making the buoyancy, an element 
of its dual, into a function. 



6 Hamiltonian formulation of Euler- Yang- Mills 

Once the Lagrangian formulation of a theory is known, one usually passes to the Hamil- 
tonian formulation by a Legendre transformation, if the Lagrangian function is non- 
degenerate. Unfortunately, in our case, this is not possible because the Legendre trans- 
formation is not invertible, as we have already seen when studying the Maxwell equa- 
tions. The trouble is that the Lagrangian function does not depend on To deal 
with this, we shall work with a new Lagrangian function obtained by eliminating ^0 
from (15. Sp . For this new Lagrangian function, the Legendre transformation is invertible 
and we can deduce the associated Hamiltonian formulation. However, in this process, an 
equation gets lost, namely. Gauss' Law 5-^£ = — Q in fl5.17p . This equation will be recov- 
ered as a conservation law of the momentum map associated to the gauge transformation 
group. We begin by quickly recalling some facts about the Hamiltonian semidirect prod- 
uct reduction theory. 



6.1 Hamiltonian semidirect product reduction with parameter 

Let S := G (§)V he the semidirect product defined at the beginning of ^ The lift of right 
translation of S on T*S induces a right action on T*G x V* . Let Q be another manifold 
(without any G or V^-action). Consider a Hamiltonian function H : T*G x T*Q xV* ^M. 
right invariant under the S'-action on T*G x T*Q x V*\ recall that the S'-action on 
T*Q is trivial. In particular, the function := -f^|T*Gxr*Qx{ao} • T*G x T*Q — >• 
M is invariant under the induced action of the isotropy subgroup GaQ := {g & G \ 
p*ao = flo} for any oq G V* . The follow ing theorem is an easy consequence o f the 
semidirect product reduction theore m (see Marsden. Ratiu. and Weinstein |l984| ) and 



the reduction by stages method (see iMarsden et all |2nn7j ). 



Theorem 6.1 For a{t) e T*^t)^ and pit) := T* Rg(^){a{t)) G Q* , the following are 
equivalent: 

i («()!:), g(t),p(t)) satisfies Hamilton's equations for on T*(G x Q). 

ii The following system of Lie-Poisson equations with parameter coupled with Hamil- 
ton's equations holds on s* x T*Q: 

-^(/i, a) = -ad/«,, 6h\ip,a) = - ( ad*sh p + ^ oa, ) , a(0) = 
at {jji^j^) \ 4m da dp J 
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and 

dq^ dh dpi dh 
dt dpi^ dt 

where s is the semidirect product Lie algebra s — Q®V. The associated Poisson 
bracket is the sum of the Lie-Poisson bracket on the Lie algebra s* and the canonical 
bracket on the cotangent bundle T*Q, that is, 

{/ 9}ipi aqp)=(fi ^1 \ + /a - + - 

\ ' _b\x b\x_ I \ ^SaS/j, 6a6iJ,/ dq^ dpi dq^ dpi 

For example, one can start with a Lagrangian as in the previous section, sup- 
pose that the Legendre transformation FLa^ is invcrtible, and form the corresponding 
Hamiltonian = E^q o FL~^^, where is the energy of L^g. Then the function 
H : T*G X T*Q x V* ^ IR so defined is -S-invariant and one can apply this theo- 
rem. This is the method we shall use below to find the Hamiltonian formulation of the 
Euler- Yang-Mills equations. 



6.2 The Hamiltonian 

Recall that we identify the cotangent space T^Aut{P) with the space of G-invariant 
one-forms on P along t/j e Aut{P). The duahty is given by 

{M^,U^) := / M^U^)!!, 
Jm 

where G T^Aut{P) and U,^ G T^Aut{P). Recall that the integrand defines a 
function on M (it is independent on the fiber variables in the bundle n : P ^ M). 

For U^,V^ G T^Aut{P), the expression Kji{p) {U^{p),V^{p)) depends only on the 
class X — 7i{p). Thus K_^{U^, V^), which is a smooth function on P, does not depend on 
the fibers and hence induces a smooth function on M. Therefore the integral 

Jm 

is well-defined. Moreover, the definition of Kj( immediately implies the equality 

/ K^iU^,V^)i,= [ gm],[V^])i,+ [ ^(AiU^),AiV^))i,. 
Jm Jm Jm ^ ' 

Similarly, for M^,N^ G T^Aut{P) we can define the integral 

Jm 

where denotes the dual metric induced on T*P by the Kaluza-Klein metric. 

The Hamiltonian for the Euler- Yang-Mills equations is defined on the cotangent 
bundle T*{Aut{P) x Conn{P)) and is given, for (p, s) G J^{M) x J^{M), by 

H^,,s){M^,A,y)^l [ -K*^{M^,Mi,)i2+ [ pe{p{J^r\s)i, (6.1) 
^ Jm P Jm 

+ 1 I WEfp+l I WBfp. 

^ J M ^ J M 
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This Hamiltonian is obtained by Legendre transforming the Lagrangian i^(p,s) in the case 
the variable is absent. Indeed, we have 



M4p) := ¥L{U^){j>) = p{x)Ka{p){U4p), ■), x = n{p). 

By Theorem 16. Hamilton's equations for i^(p,s) are equivalent to the Lie-Poisson 
equations on the dual of the semidirect product Lie algebra aut(P) (s) (^(M) x J^{M)), 
together with the standard Hamilton equations on T*Conn{P), relative to the reduced 
Hamiltonian h given on (out(P) (S) (J^(M) X J^(M)))* X T*Conn{P) by 

+ U WEfp+U II5IIV 

By the Legendre transformation [/ M = Kj^^{pU, ■), these equations are equivalent to 
the equations (15.171) with = but without the nonabelian Gauss equation 

S^£ = -Q. 



6.3 The momentum map of the gauge group 

This last equation is obtained by invariance of the Hamiltonian under gauge transfor- 
mations. Indeed, consider the action of the gauge group given for t] G Qau{P), by 

(tp^A) ^ {r]-^ o'ilj,7]*A). (6.2) 

The cotangent-lift of this action leaves the Hamiltonian invariant. So, the associated 
momentum map, which is computed in the following lemma, is a conserved quantity. 

Lemma 6.2 The momentum map associated to the cotangent lift of the gauge group 
action is 

J(M,^„ A,y) = cj (^6^y - A {[j^'') o ^-i)) G 0au(P)* ~ 0au(P), 

where G T^Aut{P) is such that = Kj\^{y^,-) and a : JFg(P, g) gau(P) is 
defined in (12.81) . 

Proof. We will apply the formula J(ag)(^) = {aq,^Q{q)), which gives the momentum 
mapping associated to a cotangent lifted action of a Lie group G on a cotangent bundle 
T*Q. In our case we have G = Qau{P), Q = Aut{P) x Conn[P) and for ^ = a{f) G 
0au(P), the infinitesimal generator is given by (see (13.101) ) 



exp(tO"'o^,exp(tO*^) 
(-eoV^,d^/) G r(^,^)(^izt^(P) X Conn{P)). 



t=o 
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Thus, changing variables in the third equahty below, using (12 ■4p . we get 

(j(M^, A, y),i) = ((M^, A, y), (-e o ta, a, d-4/)) 

/ {91) (y.^f)y^ 



M 



M 



K 



A 



M 



J^'^)V^o^lJ-\i]^^+ I 7(W,/)/i 



M 



Since [^] = 0, the first term can be written as 

^(A((j^~')v^o^-A,A{0\f^. 



M 



Thus, using the pairing (I2.1ip and the identity Ao a = idy^^^p^g), we get 



{3{M^,A,y),0 



M 



A/ 



+ / 7 5-^3^,^(0 U 



a(6-^y-A(^(^J^'')v^o^'^^^,^ 



When = FL(f/^) = Kj({pqU^, ■), is a solution of Hamilton's equations associated 
to the conservation law J(M^,^, 3^) = gives 



The definition of the charge density Q (see (15.131) without ^o), the identities U^oip ^ 



U, l^Jip j (po ° ) = Pi ^iid the notation 8 = —y, gives 

Q = -5^£. 

The following theorem summarizes the results of the present section. 

Theorem 6.3 Let (M^,^, 3^) be a curve in T*{Aut{P) x Conn{P)) and consider the 
induced curve (M,^, 3^) G aut(P)* x T*Conn{P) given byM:= {Jip) o ip-'^. Then 
(M^,^, 3^) is a solution of Hamilton's equations associated to the Hamiltonian if(p(,^so) 
given in (16.11) if and only if (M, A, y) is a solution of the system 

f dv 1 « 1 

_ + V.i; = -7 (g, E{-) + 5(-, v)Y - - gradp, 

dp ds 

^ + div(/w) = 0, p(0) =po, — + dsiv) = 0, s(0) = So, 



dt 

^ + V^Q + Qdivv = 0, 



dt 



(6.3) 



d£_ 

'at 

dB 
'dt 
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where we use the same notations as in Theorem \5.4\ except that here v and Q are given 
in terms o/M by 



M 

V = [U] and Q = A{pU) where U = K*^(—, 



(6.4) 



Conservation of the momentum map associated to the gauge transformations gives the 
equation 

One can adapt this theorem to the incompressible and homogeneous case. 



6.4 The Poisson bracket 

From Theorem 16. II we know that the Euler- Yang-Mills equations (16. 3p can be written as 



with respect to the Poisson bracket on [aui{P)®{T{M) x T{M))]* x T*Conn{P) 



{f,g}{m,p,s,Ay) 
+ 



M 



M 



I A 


\5f 






_6M 






'5g] 


-< 




_6M_ 





(6.5) 



Sf_ 
6M 



(5f 


\5g] 




\5f] 


)) 


\6s 


_6M_ 


_6M_ 





M 



6A 6y^ 



We can obtain this bracket and the associated Harn i lton e quations (16.31) alternatively 
by a reduction by stages process (see iMarsden et.~aD \2mi\ ). The symplectic reduced 
spaces are of the form O x T*Conn{P), where (9 is a coadjoint orbit of the semidirect 
product S := Aut{P) (S) (J^(M) x J^(M)). 

If the principal bundle is trivial, the automorphism group is the semidirect prod- 
uct V{M) (§)J^{M,G) of the diffeomorphism group of M with the group of G- valued 
functions on M. In this case the first term can be written more explicitly by taking 
advantage of the internal structure of Aut(P \ and we recover (up to sign conve i itions ) 
the Poisson bracket given in equation (38) in Gibbons. Holm, and Kupershmidt [l98^. 



6.5 Summary 

We comment now on the structure of the equations of motion (16.31) and the Poisson 
bracket ( 16. 5p . Note that in ( 16. 3p there is an evolution equation for the gauge charge 
Q but that the functions for which the Poisson bracket (16. 5p is defined seem not to 
depend on Q. The explanation of this fact is given in Theorem 16. 3t the discussion below 
summarizes briefly the key results and comments on the structure of both the equations 
and the Poisson bracket. 
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(1) The equations for v and Q are the "components" of a single equation: the Lie- 
Poisson equation on the dual of the Lie algebra of the automorphism group. The 
true variable is the fluid momentum M which defines both the Eulerian velocity 
V and the gauge-charge Q by using (16 .4^ . Conversely, given p, t>, and Q, the fluid 
momentum M is found by putting M := pK_^{U, ■) = pg{v,Tn{-)) + 'y{Q, A{-)), 
where U := Hor^oti + ^cr(Q) (recall that for any U G aut(P) and A G Conn{P) 
we have the identity U = Hor_4o[f/] + a{A{U))). In other words, the Lie-Poisson 
equation for M is equivalent to two equations: the equation for v and the equa- 
tion for Q. This is the usual Kaluza-Klein point of view, namely, Q and v are 
constructed from M and vice-versa. 

(2) The Q-equation looks like advection equation. To see this, recall t hat Q G P/Ad P ) 
and that 6^ : fii(M,AdP) P(AdP) (see Definition 4.2.8 in Isieeckei] |l98l| l 
Defining 

div^(Qi;) := -d-^iQ (g) v^), v G X(M), 

where the one form G f2^(M, AdP) is given by {Q®v'"){ux) := Q{x)gx{v{x),Ux 

for any Ux G T^M, one easily deduces the formula 

dw^{Qv) = V^Q + Qdwv 

which allows us to write the Q-equation in the form 

Q + div'^(g^;) = 0. 

However Q is not advected since its evolution is not given by the pull back of the 
flow of the velocity field v. Note that in this equation A is itself a variable that is 
time dependent. 

(3) The equations for p and s are usual advection equations for a density and a function 
that appear in the Lie-Poisson equations of a semidirect product. 

(4) The equations for S and B are Hamilton's equations for the conjugate variables 
iA,y) G T*Conn{P). 

(5) The equation = —Q comes from momentum conservation associated to gauge 
group symmetry and d'^B = is the Bianchi identity for the connection A and its 
curvature B. 

(6) The Poisson bracket (16. 5p contains two types of terms: the first three are Lie- 
Poisson for a semidirect product and the fourth is the usual bracket on T*Conn{P). 
However, note that the first summand in (16. 5p gives rise to an evolution equation 
for M which, as we discussed above, is equivalent to two evolution equations, 
one for v and and another one for Q. If the bundle is trivial, one can make the 
formulas (16. 4p more explicit, as we shall see below when we carry out one more 
reduction. Note also that the Poisson bracket (16.51) is a product bracket; there 
is no coupling between the semidirect product fluid variables (M, p, s) and the 
Yang-Mills field variables (.4., 3^). The coupling in the equations is exclusively due 
to the Hamiltonian (16. ip . 
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6.6 The second reduction 



Note that right translation in the group 5* = Aut{P) (§) {J^{M) x ^(M)) on itself and 
the action of Qau{P) on Aut{P) x Conn{P) given by f l6.2p commute if one views them 
as actions on S x Conn(P). There fore, by the general theory of commuting reduction by 
stages (see iMarsden et. al. since the momentum map associated to the gauge 



group action is ^Mt(P)-invariant, it induces a momentum map J^* on s* x T*Conn{P) 
which restricts to a momentum map Jq on the reduced space O x T*Conn{P). Here 
s := aut(P) (S) (^(M) x JF(M)). A direct computation shows that the momentum map 
Js* : 5* X T*Conn{P) — > gou(P)* has the expression 

J,* (M, p,s,A,y) = cr - A{V)) , (6.6) 

where V G aut(P) is such that M = K^{V, ■), and a denotes the map defined in fl2.8p . 
The gauge group action induced on s* x T*Conn{P) and O x T*Conn[P) is given by 

(M, p, s, A, y) ^ (Ad; M, p, s, r]*A, r]*y). (6.7) 

Using the notation S := (M,p, s) G O, it can be written as 

{S,A,y)^{Adl^^o,o)S,V*-A,v*y). (6.8) 

This action is simply the diagonal action given on the first factor by the coadjoint action 
of the subgroup Qau{P) of S = Aut{P) (§) {J^{M) x JF(M)), and on the second factor by 
gauge transformations. Note that when the center Z{G) of the group G is trivial, then 
the transformation A t— > ri*A is free. In this case, the reduced action (16.81) is also free 
and the second reduced symplectic spaces 

J^i(N)/to(P)N, Negau{Py, 

have no singularities. 

By the reduction by stages process, the reduced spaces 3Q^{N)/Qau{P)j<! are sym- 
plectically diffeomorphic to the reduced spaces obtained by a one step reduction from 
the cotangent bundle 

T*{S X Conn{P)) 

with respect to the product of the two cotangent-lifted actions. Note that these reduced 
spaces are, up to connected components, the symplectic leaves in the Poisson mani- 
fold J-}(N)/Gau{P)T^. This is a straightforward consequence of Theorem lO.l.l(iv) in 



Ortega and Ratiij |2004| . because the optimally reduced spaces are, up to connected 
components, precisely the symplectically reduced spaces for every leaf. 

Note that the Euler- Yang- Mills equation (16.31) projects to the reduced space at zero 
momentum 

3o\0)/GauiP). (6.9) 

The general case corresponds to the Yang-Mills fluid with an external charge N. 

In order to obtain the reduced Poisson structure concretely, we will identify the space 
s* X T*Conn{P) with a space on which the gauge action is simpler. This identification 
is given in the following proposition. 
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Proposition 6.4 Consider the group K = V{M) (§) {J^{M) x J'(M)) and denote by t* 
the dual of its Lie algebra. There is a gauge-equivariant diff'eomorphism 

i:s*x T*Conn{P) ^ T x TaiP. 0*) x T*Conn{P), (6.10) 

given by 

2(M, p, s, A, y) := ((Hor^)* o M, p, s, JJ o M, A, -y) =: (n, p, s, z/, A, £), 

where the gauge group acts on s* x T*Conn{P) by the action (16.81) and on t* x T*Conn{P) 
only on the factor J^c{P,d*) x T*Conn{P) by the right action 

(z/, A, £)^{iyo 7], r]*A, r]*£). (6.11) 

Moreover, the image of the level set J~*^(N) by the diffeomorphism i is 

{(n, p, s, u,A,S)\iy + ^{6^8 + /,-) = 0}, 

where N G Qau{P) and f G J^g{P-, q) is such that o{f) = N. Thus J~,"^(N) is diffeomor- 
phic to r X T*Conn{P). 

The map J : T*P — > 0* denotes the momentum map J(ap)(^) := {otp,^p{p)), and 
(Hor^)* denotes the dual map of the horizontal-lift Hor_4 : TM — >■ TP with respect to A. 

Proof. We first prove that i is injective. Suppose that (M, p, s, A, y), (M', p', s\ A', y') G 
s* xT*Conn{P) have the same image under i. We clearly have (p, s, A, y) = (p', s', A', y'). 
Therefore we have (Horyi)* o M = (Hor^)* o M' and J o M = J o M'. This imphes that 
M = M'. The map i is clearly surjective and hence invertible, its inverse being given by 

i^'^{n, p, s, u, A, 8) = (7r*n + A*^, p, s, A, —8), 

where A* : Q* T*P denotes the dual map of A. It follows that i is a diffeomorphism. 
To prove gauge-equivariance, it suffices to show that for all rj G Qau{P), 

(Hor^M)* o = (Hor^)* o M and Jo ri*M = (J o M) o r/. 

This is a direct computation using the formulas 

T;7r(Hor^);(M(p)) = M(p) - ^(p)*(J(M(p))) 

ri*A = Adff-i oA + TLff-i o Tr] 

rj*Mip) = M(p) +T3T*%p)-.(J(M(p)))), 

where rf G J^g{P, G) is such that rj{j)) = $f^(p)(p). 

Recall that Js*(M, p, s, A, y) = a^d-^y - A{V)), where V G aut(P) is such that M = 
K_AiV, ■). So for N = a{f) the condition J,*{M, p, s, ^, 3;) = N reads S-^y - AiV) = f. 
Using that z/ = J o M = 'y{A{V), ■) and 8 = —y, we get the condition 

z. + 7(5-^^ + /,-) = 0. ■ 
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This proposition shows that the reduced spaces J~}(N)/Qau{P)]ss can be identified 
with the quotient t* x [{J^g{P,Q*) x T*Conn{P)) /Qau{P)i^}] via the diffeomorphism in- 
duced by i and given by 

[(M, p, s, z/, A, y)] ^ ((Hor^)* o M, p, s, [z/, A, S]), (6.12) 

where [■] denote the corresponding equivalence classes. 

We now compute the Poisson structure { , }' induced by i on t* x J^g{P,Q*) x 
T*Conn{P). For f,g E T{t* x J^g{P,Q*) x T*Conn{P)) we have the formulas 



^ - Hor^ o-^ + J* o 

on ou 



6M 



{91) 



5M 



n:=(Hor^)*oM, 

on 



5A ' 



C 



{91) 
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C] -uiCi Rot A 



(£) 



1/ := J oM. 



Using the equality (15. 8p . we obtain 

{/, 9yi^, p, s, i^, A £) ■= {f °h9° p, s, A, y) 



(6.13) 



n 



M 



+ 



M 



SI 6g_ 
5n' 6n 

5f 



M 



5f \ 6g 



s ( div ( -r^^ 1 — div ( -r^^ 



Ss 6n 



Ss 6n 



5p J Sn 



S9\ 5f 



Sp J Sn 



A/ 



+ 
+ 
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6A' 6S 
SI 5g_ 
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A/ 
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6£ 



5f\ 5f 



6n 



\ TF \ Hor^ o— — — Hor^ o— 



6£ 
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TT ^9 

Hor^o — 
on 



+B [ Rota o|^, Hor^ 

on on 



6n 



p. 



Hor^o — 
on 



Note that the first three terms in (16.131) represent the Lie- Poisson bracket on t*, the 
fourth and fifth terms represent the canonical bracket on T*Conn{P), the sixth term 
is the Lie-Poisson bracket on J-'g{P,Q*), and the last term provides the coupling of the 
fiuid variables to the Yang-Mills fields. 

By the general process of Poisson (point) reduction, the reduced spaces 
J;}{N) /GauiP)^ ~ r X [(J^GiP^d*) X T*Conn{P)) /GauiP)^] 
inherit a Poisson bracket { , }n given by 

{/n, ^N}N(n, p, s, [u, A, S]) := {/, ^}'(n, p, s, u, A, y), (6.14) 
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where f,g are any ^aM(P)-invariant extensions of the functions /n o ttn, (/n o ttn : 
J-*^(N) M, relative to the projection ttn : J7*^(N) J-^(N)/^aM(P)N. 

There are no exphcit formulas for the equations of motion on the Poisson point 
reduced space jQ^{0)/Qau{P) because there is no concrete realization of this quotient, 
to our knowledge. However, there is an important particular case where this is possible 
that we study next. 



6.7 The case of a trivial bundle 

We end this section by examining the case of a trivial principal bundle P = M x G 
and, more precisely, the case of the Euler-Maxwell equations which are obtained by 
ta king by G = S\ We then compare our results to those obtained for Euler-Maxwell 



m 



Marsden et.~aD |l983| . Recall that we have out(P) = X(M) (s) J'(M, g), so we obtain 



aut(P)* = X(M)*(s)J^(M,g)* = n^{M)(§)T{M,g*). For (m, z/) G aut(P)*, the gauge 
transformation fl6.7p is given by 

(m, u,p,s,A,Y) ^^ 

(m + T*ri o T*Rff~i o z/, Adi oz/, p, s, Ad.;^-i oA + TLff-i o Trf, Ad^j-i oY), 

where A := A,Y := y E Q{M, g). The equivariant diffeomorphism (16.101) is 
i(m, z/, p, s. A, r) = (m - AV, p, s, u, A, -F) =: (n, p, s, u, A, E), 
and the gauge transformation (16.111) is 

(n, p, s, z/, A, (n,p, s, Adioz/,Adj^-i oA + TL^-i oTf/,Ad^-i oE). (6.15) 

Recall that the relation between the charge density Q and the momentum M is 

Q = AiV), 

where V E aut(P) is such that M = K_a{V, ■); see (16.41) . When the bundle is trivial, this 
relation reads 

z/ = 7(Q,-). (6.16) 
In the case of Euler-Maxwell, since G = S^, the gauge transformation is simply 

(m, z/, p, s. A, F) t-^ (m + udr], z/, p,s,A + dr], F), 

and the relation (I6.16P reads v = Q. Recall that we can write Q = p^, where q G !F{M) 
is the charge, see C oroUarv 15.51 Th i s gau ge transformation coincides with the one given 
in equation (36) in lMarsden et. al.l 1983l |. where the notation a := ^ is used. The zero 
level set of the momentum map is 



J,-.i(0) = |(m,p^,p,s,A>^) 5Y = p^,]. 



The bijection i reads 



( 1 \ ( Q Q 

i m, p— , p,s,A,Y] = m - Ap—,p, s, p— , A, -Y 
\ m /V mm 
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and the image of J5"*^(0) is 

{(n, p, s, z/, A, E) I div E = u}, 

where the notation E := i?" G X(M) is used. The gauge transformation (16.111) is 

(n, p, s, z/, A, E) ^ (n, p, s, z/, A + drj, E). 

Through the diffeomorphism i, the projection ttq : J7*^(0) 3^}{0)/Qau{P) is given by 

(n, p, s, z/, A, E) ^ (n, p, s, [A],E), 

where [A] denotes the cohomology class of the one- form A. Assuming that the first 
and second cohomology groups of M are trivial, H^{M) = H'^{M) = {0}, we get the 
isomorphism 

[A]^ B ■.= dAenl{M), (6.17) 

where r2^;(M) denotes the space closed two-forms. Thus i induces a diffeomorphism 
between J-}{0)/gau{P) and the space t* x fi^^(M) x n^{M) given by 

[m, z/, p, s, A, Y] (m — Au, p, s, dA, —Y) =: (n, p, s, B, E). 

where fi^;(M) denotes the close d 2-forms on M . This identification coincides with the 
one given in Proposition 10.1 in iMarsden et. al.l 1983 . 



Using the definition (16.141) and the bracket (16.131) . the reduced Poisson bracket on 
r X Vtli{M) X fii(M) is 



{/,^}o(n,p,s,5,E) 
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M 



SI 5g_ 
(5n' 6n 
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M 



m V 8E \ (5n 



5E \6n 



5n 5n 



p, 



and the Euler-Maxwell equations can be written as 

/ = {/,Mo, 

relative to the induced Hamiltonian h given by 

h{n,p,s,B,E) =]- I (-g{n,n) + -{5Ef]p+l pe{p,s)p 

+ 1 [ {\\Er + \\B\n P- 
^ Jm 
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Note that the function 



C(n,p,s,S,E) = ^ / -{6Ef 

^ J M P 

is a Casimir function, so an equivalent Hamiltonian is 



h{n, p, s,B,E) 



M P 



-g{n,n)fi + 



pe{p, s)p + - 

M ^ 



M 



When M is three dimensional, we can use the notations B 
Therefore the two last terms can be written as 



\Ef+\\Bf) p. 

-- and E := EK 



M 



+ 



p 



M 



,5/ 5g 



p 



P— 

M m 



5g_ 51 



9 



6f 6 



' 6n 6n 



p. 



This bracket coincides with the one derived in iMarsden et. al.l |1983l | by a direct 
computation. Note that the first line in the formu la above is the Pau l i-Born -Infeld 
Poisson bracket for the Maxwell equations (see, e.g. iMarsden and Ratiul 199J], §1.6). 



The Hamiltonian h is very simple: it is the sum of the total energy of the fluid plus the 
energy of the electromagnetic field. 

Remark. In the Euler-Maxwell case, the correspondence fl6.17p is a bijective map if 
H\M) = H^{M) = {0}. Indeed, for B,B' such that dB = dB' = 0, we have B = dA 
and B' = dA', therefore if B = B' we have A = A' + dr], that is, [A] = [A']. 

This fact does not generalize to the case of a nonabelian principal bundle, trivial 
or not: there exist gauge inequivale nt connections (even on x G) with t he same 
curvatures and holonomy groups; see iMontgomeryl jl986l | , iGu and Yangl |l977l | , iMostow 
|l98nj . 



7 The Kelvin-Noether Theorem 

The Kelvin-Noether theorem is a version of the Noether theorem that holds for solutions 
of the Euler-Poincare equations. An application of this theorem to the ideal compressible 
adiabatic fluid (see (II. 2p ) gives the Kelvin circulation theorem 




where 7j C M is a closed curve which moves with the fluid velocity v, and T = de/ds is 
the temperature. 



7.1 Kelvin-Noether Theorem for semidirect products 

In order to apply this theorem to the Yang-Mills flu id, we recall some facts about the 



Kelvin-Noether theorem for semidirect products (see iHolm. Marsden. and Ratiul |1998 
for details). 
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We start with a Lagrangian Lag depending on a parameter oq G V*, as at the 
beginning of ^ We introduce a manifold C on which G acts on the left and suppose we 
have an equivariant map JC : C x V* ^ q** , that is, for all g E G,a E V*,c E C, we have 



Theorem 7.1 (Kelvin-Noether.) Fixing cq G C, let ^{t),q{t),q(t),a(t) satisfy the 
Euler-Poincare equations and define g{t) to he the solution of g(t) = TRg(^t)(,(t) cind, 
say, g{0) = e. Let c{t) = g{t)cQ and I{t) := I{c{t),^{t),q{t),q{t),a(t)). Then 



7.2 The Kelvin-Noether Theorem for Yang-Mills fluids 

In the case of the Yang-Mills fluid, we shall choose for the abstract Lie group G above, 
the automorphism group Aut{P) and we let C = {c G J-'{S^, P) | vr o c G Emb(S'^, M)}, 
where Emb(S'^,M) denotes the manifold of all embeddings of the circle in M. The 
left action of Aut{P) on C is given hj c ^ ip o c. The map /C is defined by 



{IC{gc,p*g{a)),fi) = (/C(c,a), Ad*/i), 



where gc denotes the action of G on C. 

Define the Kelvin-Noether quantity / : C x g x TQ x V* ^ R hj 

I{c, ^, q, g, a) := ( /C(c, a), ^(^, g, g, a) ) . 




— oa 
oa 
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A change of variables in the integral shows that /C is equivariant, that is. 



o c, {J^~\p o ^-1), s o ^-1), M) = (/C(c, (p, s)), Ad; M). 



Using the Lagrangian / given in (15.61) . we have 



61 
6U 



{U,Ao,Ao,A, A, (p, s)) 



^*iP9i[U],-))+l{Q,A{-)). 



Therefore the Kelvin-Noether quantity is 



I{c,U,Ao,Ao,A, A, (p,s)) 




i 



1 



■7(Q,^(-)), 



where c := vr o c G Emb(5'^, M). 

On the other hand, using the equality 
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we get 



/C(c, (p,s)), 
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cds 



Thus, by Theorem 17.11 the Kelvin Circulation Theorem for the Yang-Mills fluid is 

1 



d 
dt 



L"' It 



It 



Tds, 



(7.1) 



It 



where 7t := (^9^ o cq C P, 7^ = o cq C M is a closed curve which moves with the fluid 
velocity v := [U], and T := de/ds is the temperature. 

When the principal bundle is trivial, formula (17.11) reads 



d 
dt 



Tds. 



For the Euler-Maxwell fluid consisting of particles of mass m and charge q, since Q = p—, 
the second integral vanishes and we get 



d_ 

dt 



It 



m 



Tds, 



It 



which coincides with formula (7.37) in lHolm. Marsden. and Ratiul 1998 
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